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1 .  INTRODUCTION 


This  is  an  interim  technical  report  in  which  is  outlined  the 
accomplishments  obtained  to  date  by  the  Department  of  Electrical 
Engineering,  University  of  Mississippi,  in  an  investigation  of  the 
electromagnetic  properties  of  the  so-called  ground  stake  antenna 
undertaken  for  the  Naval  Ocean  Systems  Center  under  contract 
N66001-82-C-0045.  The  report  describes  progress  made  to  date  in 
three  facets  of  the  overall  research  project:  (-L)  formulation  and 
numerical  analysis  of  integral  equations  for  the  problem  of  determining 
the  characteristics  of  a  coax-fed  monopole  antenna  in  a  lossy  medium, 

(■LI)  formulation  of  the  Integral  equation  for  a  cylindrical  antenna 
which  resides  in  two  contiguous  half  spaces,  and  (■ LLL )  development  of 
the  measurement  apparatus  and  techniques  for  the  experimental 
investigation  of  the  properties  of  a  cylindrical  antenna  which  resides 
in  two  contiguous  half  spaces. 

In  Section  2  is  presented  a  detail  discussion  of  the  derivation 

of  the  coupled  integral  equations  for  a  monopole  in  a  lossy  medium 
% 

that  is  fed  by  a  coaxial  waveguide  through  a  ground  plane.  The  equations 
involve  Integrals  whose  integrands  are  singular  and  infinite  series 
whose  terms  contain  Bessel  functions,  so  to  solve  these  equations 
numerically  requires  great  care.  Therefore,  in  this  section  several 
analyses  are  described  that  enable  one  to  render  all  terms  in  the 
equations  amenable  to  efficient  and  accurate  numerical  computations. 

Section  3  is  devoted  to  a  description  of  the  formulation  of  the  integral 
equation  for  the  ground  stak?  antenna.  Because  the  ground  stake 
cylindrical  antenna  Is  partially  in  one  medium  and  partially  in  another. 


Sommerfeld  integrals  are  incorporated  in  Green's  functions  which  appear 
in  the  integral  equation.  At  best,  Sommerfeld  integrals  are  difficult 
to  evaluate  efficiently  and  accurately  'so  significant  attention  is 
given  to  the  developemnt  of  equations  in  a  form  that  can  be  solved 
numerically  with  high  efficiency.  As  opposed  to  double  integration 
which  one  must  confront  in  the  Green's  functions  for  the  ground  stake 
antenna  when  one  derives  these  Green's  functions  from  the  traditional 
Sommerfeld  integrals  for  an  elementary  dipole,  it  has  been  possible 
by  taking  a  fresh  approach  to  obtain  Green's  functions  possessing  only 
single  integrals.  A  brief  study  of  this  integral  equation  reveals  that 
it  can  be  solved  numerically  rather  readily. 

In  Section  4  is  found  a  detailed  description  of  the  apparatus 
and  procedures  presently  under  development  for  experimentally  studying 
the  properties  of  the  ground  stake  antenna.  Construction  details  of 
the  tank  in  which  the  antenna  is  to  be  placed  for  measurements  are 
presented..  Also,  the  various  pieces  of  apparatus  that  have  been 
constructed  for  the  measurement  of  the  properties  of  water  are  described. 
Finally,  the  measurement  techniques  to  he  employed  are  outlined  briefly. 
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2.  FORMULATION  AND  ANALYTICAL/NUMERICAL  PROCESSING  OF  COUPLED  INTEGRAL 
EQUATIONS  FOR  COAX-FED  MONOPOLE  ABOVE  A  GROUND  PLANE 

In  this  section  is  found  a  derivation  and  an  analysis  of  the 
coupled  integral  equations  for  a  coax-fed  monopole  above  a  ground 
plane  (Fig.  1  ).  The  monopole  resides  in  a  lossy  medium  of  infinite 
extent.  These  equations  and  the  numerical  techniques  developed  for 
solving  them  provide  the  basis  upon  which  one  may  found  an  accurate 
analysis  of  an  antenna  in  a  lossy  medium. 

The  monopole  antenna  illustrated  in  Fig.  1  is  the  extension  of  the 

center  conductor  of  a  coaxial  waveguide  whose  outer  conductor  terminates 

at  the  conducting  ground  plane.  The  monopole  is  taken  to  be  a  conducting 

tube  (open-ended)  with  a  vanishingly  thin  wall.  It  is  fed  or  driven  by  a 

time-harmonic  signal  (e^Ci't)  from  the  coaxial  guide  and  receives  its 

excitation  through  the  annular  aperture  where  the  coax  opens  into  the 

ground  plane.  The  monopole,  is  of  height  h  and  radius  H  and  it  resides 

in  a  semi-inf inite,  homogeneous  medium  characterized  by  (y,e)  with 

e  (=»  SqEj.— j  allowed  to  be  complex  to  account  for  losses.  The 

radii  of  the  inner  and  outer  coax  conductors  are  a  and  b,  respectively, 

and  the  medium  in  the  annular  region  between  the  conductors  may  be 

uniform  in  which  case  it  is  characterized  by  (yc*ec)  or  ic  ™ay  be  a 

gas  (y  ,£  )  in  which  a  dielectric  (y  ,e  )  bead  of  length  L  is  inserted 
8  8  c  c 

to  provide  mechanical  support.  (Gee  Fig.  2  .)  In  the  former  case  of 
a  uniform  solid  dielectric,  L  is  made  to  approach  infinity.  As 
suggested  in  Fig.  2  ,  the  monopole  radius  may  differ  from  that  of  the 
coax  inner  conductor. 

We  wish  to  calculate  the  current  induced  in  the  monopole  and  to 
determine  the  load,  or  admittance,  which  the  monopole  presents  to  the 


Ground 


coaxial  line.  The  known  excitation  is  specified  to  be  a  TEM  wave 
traveling  in  the  positive  z  direction  in  the  coax,  whose  field 
components  are 


E 


i 

P 


9 


H . 


(1) 


where  k2  =  co2y  e  and  ri2  =  y  /e  and  where  E.  is  the  known  complex  ampli- 

tude  of  the  TEM  electric  field  incident  upon  the  annular  aperture.  If 

the  coax  is  gas-filled  and  the  center  conductor  is  supported  by  a 

dielectric  bead,  k  and  n  are  replaced  by  k  and  ri  ,  respectively.  The 

c  c  8  8 

viewpoint  taken  here  is  that  the  monopole  can  be  treated  as  a  scatterer 
which  is  illuminated  by  the  field  radiated  by  the  coaxial  aperture  in 
the  screen.  This  approach  enables  us  to  apply  aperture  and  scattering 
theory  to  the  present  analysis  and  provides  the  basis  for  the  formulation 
of  a  pair  of  coupled  integral  equations  for  the  monopole-coax  structure. 
As  an  aid  in  the  derivation  of  equations,  the  formulation  is  partitioned 
into  three  parts:  field  in  coaxial  region,  field  radiated  by  annular 
aperture,  field  scattered  by  monopole. 

Formulation  of  Integral  Equations 

Coupled  integral  equations  are  formulated  in  this  section  for  the 
unknowns  I  and  E^  where  I  is  the  total  axial  current  on  the  monopole  and 
Ep  is  the  p  component  of  the  electric  field  in  the  annular  aperture  A 
through  which  the  coaxial  guide  and  monopole  antenna  are  coupled.  These 
integral  equations  are  derived  by  requiring  (4.)  that  the  <J>  component 
of  the  H-field  in  the  coax  be  equal  to  that  in  the  exterior  region  when 
both  are  evaluated  in  the  limit  as  the  points  of  observation  approach 
the  annular  aperture  and  (-c-t)  that  the  electric  field  tangential  to  the 
monopole  and  to  the  conducting  plane  be  zero.  For  this  purpose  the 
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£ 

magnetic  field  is  expressed  in  the  coaxial  guide  as  a  function  of  Ep 
and  of  the  excitation  and  in  the  exterior  region  as  a  function  of 

E3  and  I. 

P 


Magnetic  Field  in  Coaxial  Guide 


.co 


The  magnetic  field  in  the  coaxial  guide  can  be  written  in  terms 
of  the  coax  excitation  of  (1)  and  aperture  field  E  as 


CO  sc  c 

(p.z)  -  H^c(p,z)  +  Hj(p,z) 


(2a) 


where 


H^(p,z) 


b 

» 

Ea(p')GC(p,z;p’)dp' 

a 


(2b) 


and  in  which 


H^(p,z)  *  2 


0  . 

— -  cos  k  z 
n  p  c 

C 


(3a) 


is  the  so-called  short-circuit  magnetic  field  which  would  exist  in  the 
uniform-dielectric  coax  if  the  aperture  at  z  *  0  between  the  coaxial 
guide  and  the  exterior  region  were  short-circuited.  If  the  coax  is 
gas-filled  with  a  bead-supported  center  conductor,  the  short-circuit 
magnetic  field  in  the  region  -L  <  z  <  0  would  be 
E, 


jk  L 
,  8 


„3C.  .  »  ^0  _ - _ 

"<J>  p  *  n  cos  k  L  +  j  n  Sin  k  L 

g  c  c  c 


cos  k  z,  z  e(-L,0) .  (3b) 

c 


In  the  case  of  the  gas-filled  coax,  (2a)  is  valid  only  in  the  bead 
material  (-L  <  z  <  0) ,  and  it  is  assumed  that  higher  order  modes  excited 
at  the  discontinuity  at  z  ■  0  are  vanishing  small  at  z  ■  -L.  The 
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Green's  function  G 


c 


is  determined  by  standard  methods  [  1  ]  to  be 


GC(p,z;p') 


where 


,  jk  z 
— C 

n  p£nB 

c  a 


-j  2k  (z+L) 
1  -  T  e  c 

-j2k  L 
1  +  r  e  c 


Z£(-L,0) 


^ 

_S _ £ 

na  +  ’V 
g  c 


(4) 


(5) 


where 


yp) 


N0(ktqa)  J0(ktqP)  *  J0(ktqa)  W>  ’ 


(6) 


and  where 


k2  <  k2 
tq  c 


(7) 


in  which  k  ,  q*l,2,...,  is  the  qtk  root  of  4>  (b)*0.  In  (6),  J  and 
tq  b  V 

N  are  the  vtk  order  Bessel  and  Neumann  functions  and  the  norm  N2  in 
V  q 

(4)  is 


The  Green's  function  Gc  of  (4)  is  valid  for  both  the  uniform-dielectric 
coax  and  the  gas-filled  coax  with  bead-supported  center  conductor;  in 
the  former  case,  f  is  set  equal  to  zero  in  the  first  term  of  (4). 
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Field  in  Exterior  Region 

The  magnetic  field  in  the  exterior  region  is  determined  as  the  sum 

H®X  =  H*  +  H*  (9) 

<}>  <+> 

where  is  that  radiated  by  the  annular  aperture  in  the  absence  of  the 
<P 

g 

monopole  and  H.  is  that  scattered  by  the  monopole.  The  field  incident 
<P 

upon  the  monopole  is  that  radiated  by  the  aperture.  As  an  aid  in 
deriving  expressions  for  the  exterior  field,  the  exterior-field 
equivalent  model  of  Fig.  3  is  introduced.  In  this  figure  one  sees 
that  the  annular  aperture  A  at  z  *  0  is  short-circuited  and  that  an 
equivalent  surface  magnetic  current  of  density  where  * 

i  ■  ft  § x  *]  or  ■  -E^  is  placed  on  the  shorted  annulus  [  2 ] .  This 
model  and  its  field  are  equivalent  to  the  original  structure  and  its 
field  in  the  exterior  (z>0)  region  [  2  ].  By  image  theory  one  arrives 
at  the  final  model  of  Fig.  3  which  is  used  as  a  guide  in  exterior 
formulations. 

r  f 

In  the  absence  of  the  monopole,  the  field  (E  ,H  )  due  to  the 
magnetic  current  is  determined  from 


Er  »  -  ~  V  x  F 


(10a) 


Hr  -  -j  A  (k2F  +  VV  •  F) 


(10b) 


where 


*  *  11 


,  -Jkl £-p* | 

2M.  (p1  )$’  - -  dS' 

♦“  -  lr-e’1 


or  where  F  ■  F .$  with 

<P£ 


9 


F  (p,z) 

4> 


£ 

27T 


*  a  f*  e-jk[22  +  R2]?i 

p'E*(p')  !  coscf)'  e 
P  J 


[z  +  R'J 


21% 


d<t)'dD'  (12a) 


in  which  k2  *  u2y£  and  n2  *  u/e  and 


R2  *  [p2  +  p'2  -  2pp'cos<t»'  ]  .  (12b) 

From  (10)  and  (12)  one  readily  finds  that  the  field  components  Er 

z 

and  Hz  are 
<P 

-jk[z2  +  R2]1* 

cos0  '  — - r—  d^'dp'  .  (13) 

r  u2  +  *2r 

and 

i  .  fb  rn  ~jk[z2  +  R2]*5 

-  J  j  P'E*(P')  J  cos*'  .  <I4) 

a  -it 

The  field  (E^.H0)  due  to  the  total  axial  current  I  induced  on  the 
monopole  can  be  determined  from 

ES  -  -j  —  (k2A  +  W  •  A)  (15a) 

and 

HS  -  i  Vx  A  (15b) 

U  -  - 

where  A  has  only  a  z-component  given  by 

h  7T  -jk[(z-z’)2  +  Rj]*5 

A,(P.«)  "  4tt  J  I(z,)  -k  J  5 -  d*'dz'  06) 

-h  -7T  [(z-z')2  +  R2]** 

X 

in  which  R2  is  R2  of  (12b)  with  p*  replaced  by  4.  As  is  consistent  with  the 
X 

fact  that  the  excitation  of  the  monopole  emanating  from  the  aperture 
is  $  independent,  the  current  on  the  monopole  and  the  resulting 
scattered  field  are  also  0  independent.  Expressions  for  the  field 
components  E*  and  H®  follow  immediately: 


E^(p.z) 


19 

2trp  3p  P 


P*  p(P*) 


11 


E*(p,z) 


h 

I(z') 

« 

-h 


_1_ 

2tt 


TT  -jk[(z-z’)2  +  (17) 

— - j -  d<t>'dz' 


H®(P,z) 


_1_  3_ 
4tt  3p 


TT  -jk[(z-z')Z  +  K2  ] 
sl 


2  i*5 


ir  [(Z-Z’)2+Rz] 


.2,% 


d<t>'dz' 


(18) 


-h 


Integral  Equations 

The  final  integral  equations  for  E^  and  I  result  from  enforcement 
of  boundary  and  transition  conditions.  The  boundary  condition  is  that 
the  total  electric  field  tangential  to  the  monopole  and  ground  plane 
surfaces  be  zero.  The  tangential  E-field  on  the  ground  plane  is 
zero  by  virtue  of  the  use  of  image  theory  in  the  construction  of 

j  g 

expressions  for  E^  and  J|  and  that  on  the  surface  of  the  tubular 
monopole  is  made  zero  by  enforcement  of  the  following  equation: 

ES(*,z)  +  Er(4,z)  -  0  ,  ze(-h,h)  .  (19) 

z  z 

The  required  transition  condition  is  that  the  magnetic  field  tangential 
to  the  annular  aperture  A  in  the  exterior  region  and  that  in  the  coax 
be  equal  when  the  points  of  observation  approach  a  common  point  in  A. 
This  condition  is  ensured  upon  enforcement  of 


lim  H^°(p,z)  »  lira  H?X(p,z)  ,  peA 
z+0  ♦  z+0  * 


or  of 


lim  H*T°( 


(20) 


Eqs.  (19)  and  (20)  must  hold  simultaneously  and,  upon  utilization  of 
(2),  (3),  and  (13)— (18) ,  they  lead  to  the  coupled  integral  equations 
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below  for  E  and  I: 


-j 


477k 


d2  2  s 

-772  +  *  I(z’)  G  (*,z;z')dz’ 


-h 


♦i 


p 

b 

h 0 . 

P * Ep (P * )  Ga(p,z;p’)dp’ 

L  -  J 

0  ,  ze(-h 


p-/i 


and 


"  |p  J  z<2,>  GS(p,0;z')dz'  +  j  ^  j  p'Ep(p’)  Ga(p,0;p,)dp' 


Ep(P')  GC(p,0;p')dp'  -  H®C(p,0),  pe(a.b) 


in  which 


IT  “jk[  (Z-Z1)2  +  R2]*5 

G5(p,z;z’)  -  -j-  f  - - r -  d<J>’ 

27r  J  rr _ »n2  I  d 2 v5 


and 


G®(p 


,z;p,>  "t  l 


11  --  '  2  .  ,,21*5 


CO80 


«  e 


-TT 


-jk(zz  +  R2  ] 
[z2  +  R2]1* 


d$’  . 


For  subsequent  convenience,  the  integral  equations  are  written 
in  operator  form: 


E®[I;z]  +  E^[E®;z]  -  0  ,  ze(-h,h) 

z  z  p 


and 


+  H^CEpjpJ  -  HjlEa;p]  -  HjC(p,0)  ,  pe(a.b)  . 


,h)  (21a) 


(21b) 


(22a) 


(22b) 


below 


(23a) 


(23b) 
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I  and  Ep(»-M^)can  be  determined  by  solving  the  coupled  integral  equations 
and  from  knowledge  of  these  two  quantities  the  radiated  field  and  the 
input  TEM  reflection  coefficient  of  the  coaxial  guide  can  be  computed. 


Numerical  Solution  of  Coupled  Integral  Equations 

The  integral  equations  for  Ep  and  I  are  far  too  complex  to  be 
solved  by  any  but  a  well-conceived  numerical  method.  Due  to  the 
presence  of  G  ,  one  must  sum  an  infinite  series  comprising  te.ms 
involving  derivatives  of  which  depends  upon  Bessel  and  Neumann 
functions  whose  arguments  are  determined  by  a  solution  of  a  transcen¬ 
dental  equation  for  each  q.  In  addition  the  following  terms  in  the 

integral  equations  involve  integrals  whose  integrands  may  be  singular: 
s  r  C 

E  (A,z),  E  (A,z),  lim  H.(p,z),  and  lim  H  ‘  (p,z).  A  technique  has 

z+0  z+0  ™ 

been  developed  for  solving  numerically  the  integral  equations  and 
schemes  to  handle  all  the  difficulties  enumerated  have  been  devised. 

This  technique  is  described  below. 

The  current  I  on  the  monopole  and  its  image  as  well  as  pE^  in 
the  annulus  are  represented  by  linear  combinations  of  pulses  II  of  the 
forms  N 

s 

I(z)  -  l  I  II  (z)  (24a) 

,  n  n 
n-1 


and 


Na 

pE*(p)  -  l  Vn  nn(p)  (24b) 

p  n«l 

where  I  and  V  are  unknown  constants  to  be  determined,  N  is  the 
n  n  s 

number  of  non-zero  pulses  in  the  interval  (-h,h),and  Nfl  is  the  number 
of  pulses  in  the  Interval  (a(b).  The  pulse  functions  are  defined  as 
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v<> 


(  1,  c  e  <C  -AC/2,  c  +  AC/2) 
n  n 


0,  otherwise 


(25) 


where  C  is  the  pulse  center  and  A£  is  its  width.  On  (-h,h)  the 
n 

ntk  pulse  center  is  at  z^  =  -h  +  nAz  and  the  pulse  width  Az  has  value 

Az  ■  2h/ (N  +1)  while  on  (a,b)  the  ntk  pulse  center  is  at  p  =  a  +  Ap(n-^) 
s  n 

and  the  pulse  width  Ap  has  value  Ap *  (b-a)/N  .  In  keeping  with  the 

3 

boundary  condition  I(±h)  =0,  "half  pulses"  are  placed  on  the  subintervals 

(-h,-h+Az/2)  and  (h-Az/2,h),  and  their  coefficients  are  set  equal  to 

zero.  Eq.  (21a)  is  tested  [3]  with  triangles  A  »  m*l,...,N  ,  defined 

ms 

by 


A(z) 

m 


z  -  z 


1  - 


Az 


otherwise 


(26) 


to  obtain 


N  N 

s  a 

y  1  z  +  y  v  r 

L.  n  nm  *•  n 
n*l  n*l 


sa 

mn 


-  0  ,  m*l ,2, . . . ,N 


(27a) 


where 


Z  -  +  j 
mn 


and 


h 

'  * 

h 

JL.  f  / 
47Tk  J 

L<z> 

in 

d2  » 

dl?  +  k 

|  nn(z')  K(z- 

-h 

-h 

h 

mm. 

b 

f 

mm 

J  Vl> 

hD 

J  nn<P’> 

Ga(p,z;p')dp’ 

-h 

a 

dz  . 

p-4 


(27b) 


(27c) 
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In  (27b)  K(z-z')  ■  GS(^i,z;z')  is  the  exact  kernel  [  4  ]  of  cylindrical 

antenna  and  scatterer  theory.  Eq.  (21b)  is  tested  with  delta  functions 

or  is  simply  enforced  at  match  point  c  located  at  the  centers  of  pulses 

tn 

to  arrive  at 


N  N 

l  I  raS  +  l  V  [Y  -Yc  ]  =  HSC 
n  mn  n  mn  mn  m 

n*l  n-1 


(28a) 


sc 

where  H 
m 


<J>  m 


ids1)2|..«i 

where  p  ■  a  +  Ap  (m-^s) ,  and  where 


r,as 

mn 


n 

||5  l  "„<*’>  <=3(p,0;Z')dzi] 
-h 


(28b) 


P  *  P 


m 


Y _  -  j 


mn 


o 

-  f  n  (pf 

n  J  n 


)  G  (p  ,0;p')dp'  , 

m 


(28c) 


and 


mn 


n  (p*)  i  gc(P  .ojpOdp' 

n  p  m 


(28d) 


In  view  of  the  properties  of  the  pulse  function  II  of  (25)  and  of  the 

S3 

triangle  function  A  of  (26),  Z  and  T  of  (27)  can  be  simplified 

m  mn  mn 


[5,6]  to 


Z  -  j 
mn 


Az/2 

47TkAz  {  {K(zmfl 
-Az/2 


Zn  “  « 


2[1  -  *s(kAz)2]  K(z  -  z  -O 
m  n 


+  K<  Vl  *  *n  -  «} 


(29a) 
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and 


Az 

p 

Ap/2  i 

psa  m  _  1 

1  mn  K 

{  u- 

I C | /Az] 

~P  P 

|  Ga(p,  C  +  z  ;  E,  +p  )d£ 

J  m  n 

-Az 

i 

- 

-Ap/2 

(29b) 


Integration  by  parts  twice  and  the  change  of  variable  z'  »z  +  £,  are 

n 

employed  to  transform  (27b)  into  (29a)  while  the  variable  changes 

z*z  +C  and  p'*p  +£  are  employed  to  transform  (27c)  into  (29b). 
m  n 

Also  involved  in  the  conversion  from  (27b)  to  (29a)  is  a  very  accurate 

approximation  [5,6],  By  making  use  of  properties  of  IT^  and  the 

variable  changes  mentioned  above,  the  expressions  of  (28)  for  Fas  , 

mn 

c 

Y  •  and  Y  can  be  converted  to 
mn  mn 

(30a) 

(30b) 


(30c) 

To  facilitate  subsequent  discussion  it  is  convenient  to  express 
(27a)  and  (28a)  as  a  single  matrix  equation  of  the  form 


Az/2 


C  '  -  £  [Ip  I  cs<o.o:*n«;>dC] 


-un  -°‘°m 

AP/2 

Y  «  j  £  Ga(p  ,0;p  +S)dS  , 
mn  n  j  m  n 

-Ap/2 


and 


Ap/2 

C-l  • 

-Ap/2 


(31) 
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Analytical/numerical  procedures  for  computing  the  matrix  elements  from 
(29)  and  (30)  are  presented  below. 


Analytical  Aids  in  the  Numerical  Evaluation  of  Matrix  Elements 

In  this  section  are  presented  analytical  procedures  which  greatly 

enhance  the  accuracy  and  efficiency  of  methods  for  numerical  evaluation 

of  the  matrix  elements  of  (31).  In  fact,  without  the  analysis  outlined 

below  it  is  unlikely  that  sufficient  accuracy  of  matrix  element  values 

could  be  attained  to  allow  acceptable  numerical  solutions.  Z  is 

mn 

omitted  from  the  discussion  below  since  procedures  for  its  evaluation 
can  be  found  elsewhere  [  4,6  ] . 

S3 

As  a  first  step  in  casting  T  into  a  form  suitable  for  computation 

mn 

we  point  out  that  by  making  use  of  the  analysis  presented  in  [  7 ]  one 

S3 

can  show  that  T  of  (29b)  can  be  reduced  to 
mn 


r 


sa 

mn 


Az 

-j  [1-UI/Az] 
-Az 


Ap/2 

3_ 

.  3C 

-Ao/2 


Ga(A, 


S+z  ;  P  +E)d£  d£ 
m  n 


Az 

[1-  U|/Az]fca(4,  C+z  ;  p  +Ap/2) 

L  111  n 

-Az 

-  Ga(A,  S+z  ;  p  -Ap/2)  I  d;  . 
m  n 


(32) 


Evaluation  of  (32)  is  essentially  the  same  chore  as  evaluation  of 
matrix  elements  in  cylindrical  antenna/scatterer  theory  and  it  can  be 
achieved  by  methods  developed  for  that  purpose.  In  fact  for  n*l, 

3  3 

G  (A.,  C+2  *  P_“Ap/2)  *»  G  (A,  C+z  ;  a)  which  for  A»a  is  the  exact  kernel  for  the 
m  n  m 

circular  cylinder,  while  for  n  > 1  GS  is  bounded.  Hence,  from  the 
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simple  expression  of  (32),  one  can  avail  himself  of  any  of  the  numerous 
procedures  [  4  ,  6  ]  developed  for  cylindrical  antenna/scatterer  theory 

S3 

to  evaluate  T 

mn 

Q  S 

A  form  of  T  suitable  for  numerical  computation  is  developed 
mn 

by  noting  that  an  interchange  of  differentiation  and  integration 
in  (30a)  is  valid  for  p^4.  Since  the  bracketed  expression  in  (30a) 

3  S 

is  evaluated  only  at  points  p  =  p  >  4+%Ap,  T  can  be  rewritten 

m —  mn 

as 


h  ir 


,as 

mn 


1 

8it2 


-jkD 


m 


n  (z*) (1 +  jkD  ) (p  -  4coscJ>') 
n  mm 


—  d^'dz' 


-h  -it 


m 


h  ir 

n  (z')  I  d<J>’dz' 
n  as 

-h  -IT 

where 

D  =  [p2  +  A.2  -  2p  4cos  (J)'  +  z'2]'*. 
mm  m 


(33) 


The  integrand  in  (33)  is  bounded  and  one  can  evaluate  the  integral 
numerically.  Because  of  the  small  size  of  the  coax  aperture, 
however,  the  observation  point  may  be  quite  near  the  cylinder  surface 
relative  to  the  length  of  a  source  segment  on  the  cylinder.  Thus 
the  integrand  of  (33)  may  be  very  highly  peaked  for  the  cylinder 
source  segment  in  which  |z'|  <Az/2,  To  facilitate  evaluation  of 
(33)  in  this  case  the  integral  is  expressed  in  a  form  more 
appropriate  for  numerical  calculations.  We  note  that 


*as  kD  -*6 
m 


DJ 

m 


(p  -  4  cos  4>')  = 

m 


(34) 
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3  s 

Thus  when  n  =  (N  +l)/2,  T  is  expressed  as 
s  mn 


-Az  -TT 
2 


(I  -  I'  )d4>'dz'  +-3^7  I’  d4»'dz' 

as  as  8tt2  J  J  as 

-AZ  -TT 
2 


The  highly  peaked  behavior  of  the  integrand  in  (33)  is  then  localized 

to  the  integrand  of  the  last  integral  in  (35)-  Furthermore,  integration 

with  respect  to  the  variable  z'  can  be  performed  analytically  on  I' 

as 

to  yield 


I'  dz'dcj)' 
as 


-it  ~AZ 
2 


m  n 

i  i 1  i  2 1 

=  —r-  I1  dcj>’  +  I  d<J>’ 

8it2  as  T  8tt2  as 


where 


,  Az[(p  -4)  +  4(l-cos4>')] 

^  _ m _ 

38  [  (P  -4)2+2p  4(1-0084)')  ][(p  -4) 2  +  2p  4(l-cos4> '  )  +  (~)2]h 

mm  mm  i. 


(37a) 


2 1 

I  *  k2  (p  -  4cos4> 
as  m 


' )  (fcn  +  [p2  +  42  -  2p  4cos4>'  +  (^rO2]^ 
(  2  ra  m  4 


-  -r  2n(pz  +  42  -  2p  4cos4>')  > 
z  m  ml 


(37b) 


The  integrand  I  is  slowly  varying  with  respect  to  4> '  and  can  be 
as 

•  1 f 

evaluated  numerically.  The  integrand  I  ,  however,  is  still  highly 

as 

peaked  near  4>’  ■  0.  We  therefore  subtract  the  small  argument  behavior 


of  I  under  the  integral  sign  and  add  it  to  the  expression  for 
as 

cl  S  3  s 

r  as  a  separate  integral,  so  that  T  is  finally  given  by 
mn  ran 


Az  TT 


.as 

‘ran 


1 

8tt2 


(I  -  I'  )d<J>'dz’ 
as  as 


-AZ  -7T 


A  f  (I1'  -  I1  ")d4> '  +-^7  f  I2  fd4>  ’ 
8tt2  J  as  as  8tt2  J  as 


-n 


-IT 


8n2  J 


I*"d<t>’  ,  n-(N+l)/2 

as  s 


-IT 


(38) 


where  I  ,  I'  ,  I1  ,  and  1^  are  defined  in  (33) — (37 )  and  where 
as  as  as  as 


as 


Az[  (p^-4)  +  ^4>'2] 

[ (P  -*)2  +  &)2fy(p.  -*)2  +  P  *4>’2] 
in  i.  m  m 


IT 


-7T 


Az 


Az, 


Pm[(P  „-*>  +  ^r] 

ram  l 


.  P"«+*  ,  -1 

tt  +  ■  —  tan 

in 


Tl/p  4 

in 


Pm"  K 

in 


(39a) 


(39b) 


For  n  +  (N  +  l)/2,  of  course,  T  may  be  evaluated  directly  from  (33) 
s  mn 

via  numerical  integration. 
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To  cast  Y  of  (30b)  into  a  form  amenable  to  numerical  computation, 
mn 

3L 

one  first  observes  from  (22b)  that  the  integrand  of  G  is  unbounded 

for  m*n  and  that  it  can  be  sharply  peaked  for  other  combinations  of 

m  and  n.  Thus  Y  must  be  computed  with  great  care.  To  this  end, 
mn 

we  modify  the  form  of  this  integrand  by  adding  and  subtracting 

1/R  to  obtain 
mn 


cost])*  e 


-jkR 


mn 


cos4> '  e 


-jkR 


mn 


mn 


mn 


mn 


where 


(40a) 


rL  «  P*  +  (P+O2  -  2p(p  +Ocos  <j>’  .  (40b) 

mn  m  n  m  n 

The  unboundedness  is  confined  to  the  second  term  of  (40a)  while  its 
first  is  very  slowly  varying  and  can  be  integrated  numerically  with 
ease.  Furthermore,  the  integral  over  (-n,n)  of  the  second  term  can  be 
expressed  as 


71 


4 

(p  +p  +S) 
m  n 


tt/2 


_  /1-82  sinz9 
0  mn 


d0 


(41a) 


4 

(p +P +o 

m  n 


K(B  ) 

mn 


(41b) 


where  K  is  observed  to  be  the  complete  elliptic  integral  of  the  first 
kind  [  8 ]  and  where 


6L 

mn 


4 


P  (p  H) 

m  n 

[p  +P  H)S 
m  n 


(41c) 
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(45) 


The  unbounded  part  of  the  elliptic  integral  is  added  to  and  subtracted 
from  the  integrand  of  (42).  The  integral  of  the  added  terra  becomes 
the  first  terra  of  (45)  while  the  subtracted  term  remains  in  the  integrand 
of  the  first  integral  of  (45).  Therefore,  both  integrals  of  (45)  have 
very  slowly  varying  integrands  and  can  be  integrated  numerically  with  ease. 
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The  integration  indicated  in  (30c)  can  be  performed  analytically 
c 

and,  as  a  result,  Y  becomes 

mn 


mn  „  b  p 

n  in  —  m 
c  a 


-Lfn 


K-h) 


1  -  r  e 


-J2kcL 


1  +  re 


-j2k  L 
c 


k  00 

l 

n  L 


n  L',  x 

c  q«l  q  q 


d_ 

dp 


#  (p_>  [$ 

q  m  L 


VW  -  \(P 


n  V] 


(46) 


in  which  p  ,,  *  a  +■  nAp  and  p  .  ■  a  +  Ap(n-l)  and  in  which  $  (  ) 

n-Hj  n-4  do  q  m 

“I  M 


means 


$q(p)J 


P-P 


m 


The  series  in  Y  must  be  summed  efficiently  if  the  numerical 
mn 

procedure  described  here  is  to  be  practicable.  This  is  particularly 

true  in  regard  to  YC  because  of  the  need  to  solve  a  transcendental 

mn 

equation  to  determine  the  transverse  eigenvalue  k  of  (6)  prior 

to  the  computation  of  each  series  term.  To  this  end  Rummer's 

transformation  [9]  is  used  to  accelerate  convergence.  To  employ 

this  technique,  one  determines  the  q^1  term  for  large  q  of  his 

series  and  then  finds  another  series  with  known  sum  and  the  same 

large-index  q^  term.  Even  though  Rummer's  transformation  is  applied 

to  the  series  of  the  present  problem  under  the  condition  that  the 

guide  is  below  cutoff  of  the  higher-order  modes,  k2<k2  ,  which  is 

tq 

the  case  of  major  practical  interest,  the  procedure  is  immediately 
extendable  to  the  cases  in  which  higher-order  modes  do  propagate. 

One  can  show  readily  that 


tq  q-*»  b-a 


(47a) 
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and  subsequently  that 


1  —  <$  (p  )  $  (p  )  _ ►  c  fSlnq  9mn±4  Slnq 

XqMq  dP  q(Pnr  q^Pn±%  q-*50  Cmn±4  \  q1  q3  /  ^ 


where 


^mni^  J  tt2  ? 


1  b-a 

J  .2 


b-a 


Pm  Pn±H  t/(a+Ap[m-*s])  (a+Ap[  (n-ls)*:^]) 


(48a) 


‘U.H,  ’  &  <Pm  +  pni*j  -  2*>  *  $  1  , 


and 


(48b) 


nm,tH  *  b?i  <0«  *  W  *  ^  [(«-")=«)  . 


(48c) 


In  [10]  the  series  below  and  its  closed-form  sum  are  listed: 


£  S  (a)  ■  4(a),  0<a<2ir 

q-2  q 


(49a) 


where 


~  f  s  sin  qa 

Va> '  -rf- 


(49b) 


and 


4(a)  *  -  Hn(2  sin  ^)jsina. 


0<a<2ir 


(49c) 


Notice  that,  apart  from  the  factor  C  . .  ,  S  exhibits  the  form  of  the 

ran  ±*5  q 

right  side  of  (47b)  when  q-*».  Applying  Rummer's  transformation  to  the 


series  in  Y  ,  one  arrives  at 
ran 
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1 


CD  .  \ 


mn  _  „  b  p 
H  in—  m 
c  a 


n+4 


pn-hj 


1  -  r  e 


-j2kcL 


1  +  r  e-J2k‘L 


W? 


4-  *  <p  ) fi  (p  *  (p  . )[ 

dp  i  m  J_  i  n-fV  l  n-^J 


k  °° 

—  I 

\  q-2 


1  4r  Mp_>  I*  (p -) 


XqN2  dp  W  J_  q  v n+*j 


$ 

I 


*  (P 

q 


w] 


s 

t  « 


"  Cmn+*s  Sq(0ran-^5)  "  Sq(fimn-^5) 


+  C  ,  S  (0  ,  )  -  S  (fi  ,  ) 

mn-H  q  mn-Jj  q  mn-Jj 


k  p 

c  1 


n 


-rUl4(W 


l«, 


mn+H 1 


—  C  J6 

n  ran -*5 
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which  converges  in  very  few  terms.  Absolute  values  are  used  above 
because  can  be  negative  and,  as  indicated,  (49c)  is  valid  only 


for  0<a<2ir.  Since,  for  a<0,  sin  a  -  -sin  a  ,  we  can  write 


sin  a  ■  — —  sin  la 


for  all  a 


(51) 
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This  can  be  used  to  advantage  to  rewrite  (49)  as 


?  sin  qa  a  “  ^in  q  |a|  a  x  U  , 

l  q2_i  “  —  I  — - 77  4  '  *n(2  sin  — }  slnM 

q-2  M  a  q-2  M  a 


for  -2TT<a<2tT 


S  (a) 
,  2  q 


ln(2  sin -LyL)  sin  |a|,  -2TT<a<2tr. 


(52a) 


(52b) 
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3.  FORMULATION  OF  INTEGRAL  EQUATION  FOR  GROUND  STAKE  ANTENNA 

The  analysis  of  the  behavior  of  the  current  and  charge  on  a 
straight  cylindrical  antenna  which  resides  in  two  contiguous  half 
spaces  is  a  major  portion  of  this  investigation.  As  a  model  for 
study  of  the  behavior  of  charge  and  current  everywhere — in  particular, 
at  the  wire  ends  in  lossy  media  and  at  the  interface  between  the  two 
half  spaces — the  current  and  charge  on  the  structure  shown  in  cross 
section  in  Fig.  1  is  studied.  This  section  is  devoted  to  the 
formulation  of  the  integral  equation  for  the  antenna  structure  of 
Fig.  1.  In  this  figure,  one  sees  a  circular  cylindrical  (tubular) 
antenna,  partially  in  one  half  space  and  partially  in  the  other, 
which  intersects  perpendicularly  the  planar  interface  between  the 
two  media.  In  general  the  cylindrical  antenna  does  not  meet  the 
thin-wire  conditions  and  it  is  driven  by  the  slice  generator 
illustrated.  Due  to  the  symmetry  of  the  structure  and  excitation, 
the  current  on  the  perfectly  electrically  conducting  (PEC)  antenna  is 
z-directed  and  circumferentially  independent. 

An  integral  equation  for  the  current  on  this  antenna  is 
formulated  with  care  that  it  be  amenable  to  numerical  solution. 

The  equation  incorporates  the  exact  kernel  and  special  forms  of  the 
so-called  Sommerfeld  integrals.  The  initial  phase  of  the  formulation 
is  the  derivation  of  suitable  Green’s  functions  for  the  magnetic 
vector  potential  in  the  two-media  space  of  Fig.  1  due  to  a  ring 
dipole  source.  Then  the  integral  equation — actually,  coupled 
integral  equations — is  developed  by  enforcing  the  condition  that 
the  total  electric  field  tangential  to  and  evaluated  on  the  PEC 
surface  of  the  cylindrical  tube  be  zero. 


Green's  Functions 


In  Fig.  2  is  illustrated  an  ensemble  of  elementary  dipoles 
all  on  a  circle  of  radius  p'  in  a  plane,  all  in  the  same  direction 
perpendicular  to  this  plane,  and  each  vanishingly  close  to  its  immediate 
neighbor.  This  so-called  ring  dipole  is  placed  in  a  cylindrical 
coordinate  system  with  the  center  of  the  ring  on  the  z  axis  at 
z*z'  and  with  its  plane  parallel  to  the  z»0  coordinate  plane.  The 
elementary  dipoles  of  the  ensemble  are  z-directed.  It  should  be 
clear  that  such  a  ring  dipole  source  constitutes  a  volume  current 
density  J_  given  by 

J(r)  -  z  6 (p-p ' )  6 (z-z' )  .  (1) 

It  should  be  equally  clear  that  the  magnetic  vector  potential  due  to 
J  of  (1),  with  p’  «<*  and  with  this  current  in  the  appropriately 
specified  physical  environment,  is  the  Green's  function  for  the  vector 
potential  due  to  the  current  on  the  cylinder  of  Fig.  1. 

The  Green's  function  needed  in  the  integral  equation  for  the 
cylindrical  antenna  is  the  magnetic  vector  potential  due  the  ring 
dipole  of  Fig.  2  in  the  two-media  space  of  Fig.  1.  This  vector 
potential  must  be  known  for  z  in  the  interval  (-d,h)  when  the  ring 
dipole  is  located  in  (0,h)  and  when  it  is  located  in  (-d,0).  To 
obtain  the  vector  potential  we  determine  the  solution  of  the  wave 
equation  subject  to  the  radiation  condition  with  the  ring  dipole  in 
(0,h)  and  then  with  it  in  (-d,0).  We  construct  the  solution  as  the 
sum  of  a  homogeneous  and  a  particular  solution,  and,  in  both  cases, 
we  take  the  particular  solution  to  be  the  vector  potential  due  to  the 
ring  dipole  located  in  unbounded  space,  i.e.,  homogeneous  space  of 
infinite  extent.  Hence,  in  what  follows  immediately  below,  the  wave 
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equation  is  solved  for  the  vector  potential  due  to  the  ring 
dipole  in  unbounded  space.  Next,  the  appropriate  wave  equation 
is  solved  for  the  vector  potential  due  to  the  ring  dipole  In 
region  a  and,  then,  the  appropriate  wave  equation  is  solved  for 
the  vector  potential  due  to  the  ring  dipole  in  region  b. 

The  Hankel  transform  (zero  order)  is  well  suited  as  an  aid 
in  solving  the  wave  equation  in  cylindrical  coordinates  and, 
therefore,  it  is  employed  below.  For  a  given  function  f(p)  the 
Hankel  transform  pair  is 


f(t) 


CO 

* 

pf(p)  JQ(^p)dp 

0 


f(p) 


00 

* 

Af(A)  JQ  (A.p)dA 

0 


(2a) 


(2b) 


where  f  (-t)  is  the  Hankel  transform  of  f (p) . 


Green's  Function  for  Ring  Dipole  in  Unbounded  Space.  For  the  ring 
dipole  located  in  unbounded  space  characterized  by  (y,£)  (Fig.  2), 
the  vector  potential  A  =  Az(p,z)z^  has  the  single  component  which 
satisfies 


(V£z  +  k2)  Az(p,z)  =  — y6 (p— p ' )  fi(z-z')  (3a) 

or 

4  4-  (PA  )  +  4"T  A  +  k2A  *  -  uS(p-p')  6(z-z' )  (3b) 

p  9p  z  oz^  z  z 

where  k2  ■  co2ye,  and  the  field  computed  from  Az  must  satisfy  the 
usual  radiation  condition.  Taking  the  Hankel  transform  of  the  two 
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sides  of  (3b)  (or,  equivalently,  multiplying  both  sides  by  pJg(Ap) 
and  integrating  over  (O,00))  yields 


<pAz(p-z))} 

0 


JQ(Ao)dp  + 


fa2 

3z2 


+ 


pAz(p,z)  JQ(Ap)dp 


=  -up'JqO'ip')  6  (z-z ' )  •  (4) 

In  view  of  the  properties  that  A  must  exhibit  for  p  -*■  0  and  p-*-00,  one 

z 

can  show  by  integration  by  parts  twice  that 


OO  00 

P  {p  Ip  <PAZ<P»Z))}  =  ~^\  P  VP,Z)  J0(AP)dP  (5) 

0  0 

whose  right  hand  side  is  recognized  to  be  the  Hankel  transform 

A  (A., z)  of  A  (p,z).  Substitution  of  (5)  into  (4)  leads  to  the 
z  z 

differential  equation 


dz' 


+  8 2 


Az(A,z)  =  -yp'  Jq  (Ap ' )  <5  (z-z ' ) 


(6a) 


where 


82  -  k2  -  A2 


(6b) 


The  solution  of  this  differential  equation  is 


-j 6 | z-z' 


Az(*,z)  --jf  P'J0(*P')S_ 


(7) 


from  which  A  is  found  as  the  inverse  Hankel  transform  of  A  to  be 
z  z 


Az(p,z)  -  -jy  p' 


A  J0(A.p') 


^jBlz-z' 

8 


JQ(Ap)dA  . 
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(8) 


The  form  of  the  solution  (7)  of  the  differential  equation  (6a) 
adopted  for  k^  is  such  that  a  root  6  of  (6b)  can  be  selected  in  such 
a  way  that  k^p.z)  of  (8)  represents  a  wave  traveling  away  from  the 
ring  dipole  and,  in  particular,  satisfies  the  radiation  condition  as 
1  z |  The  proper  branch  of  8  may  be  defined  by  means  of 

lm(8)  <  0  .  (9) 


Branch  cuts  which  are  consistent  with  (9)  and  which  restrict  4  to 
the  proper  sheet  in  the  4-plane  are  illustrated  in  Fig.  3.  Also 
illustrated  in  Fig.  3  is  the  integration  path  of  the  integral  of  (8). 
One  notes  that  for  k  real  the  path  must  be  deformed  into  the  first 
quadrant  to  avoid  the  singular  point  of  8  at  4  =  k,  while  in  the 
lossy-medium  case  for  which  k  is  complex  the  singular  point  is  not 
on  the  integration  path  and  deformation  is  unnecessary.  The  condition 
(9)  endows  the  integral  (8)  with  properties  which  allow  one  to  deform 
the  integration  contour  off  the  positive  real  axis  of  the  4-plane. 

If  k  is  real  and  deformation  of  the  path  of  the  integral  from  the 
real  line  is  not  to  be  employed,  8  may  be  defined  in  the  simple  way 
below: 


6  -  { 


A2  -  4 2  ,  k2  >  42 


-jA2-  k*  ,  k2  <  42 


(10) 


For  some  applications,  it  is  convenient  to  convert  (8)  to  an 
integral  over  (-00,00) .  To  do  this,  one  replaces  one  of  the  Bessel 
functions  of  (8)  by  its  equivalent  [  1  ]  in  terms  of  Hankel  functions. 
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(a) 


(b) 


Branch  points  (4. »  ±k) ,  branch  cuts,  and  integration 
path  for  (8)  in  the  -'i-plane  (a)  lossless  case,  (b) 
lossy  case. 


(11) 


J0(O  *  \  [H02)<^>  -  H^2)(Ce'jlI)]f 

replaces  the  integration  variable  by  its  negative  in  one  of  the  integrals, 
makes  use  of  Jq(~0  ■  Jq(£)»  and  combines  the  two  integrals  to  obtain 


-j6|z-z'| 

/u0(/ip*)  ^ 


P>P'  (12a) 


Az(p,z)  -  ■{ 


K  H<2)  (40*)  - 


-j8 1  z-j 


jQ(Ap)dA  ,  0<p  d2b) 


The  path  of  integration  of  the  integrals  of  (12)  and  the  branch  cuts 
(2) 

of  6  and  are  illustrated  in  Fig.  4. 

Green’s  Function  for  Ring  Dipole  in  Region  A.  In  this  subsection, 
we  determine  the  magnetic  vector  potential  due  to  a  ring  dipole  in 
one  homogeneous  half  space  which  is  in  contact  with  a  second 
homogeneous  half  space  of  different  electromagnetic  properties. 

The  two  half  spaces  are  separated  by  a  planar  interface  and  the 
ring  dipole  of  radius  p'  is  above  and  perpendicular  to  this  plane 
as  illustrated  in  Fig.  5,  where  one  observes  that  the  ring  dipole 
is  on  the  z  coordinate  axis  at  z«z'.  The  upper  half  space  is 
designated  region  a  and  the  lower  half  space  is  designated  region  b; 
quantities  peculiar  to  region  a  carry  subscript  or  superscript  "a" 
while  those  peculiar  to  region  b  carry  "b."  In  the  z  direction. 
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Fig.  4. 

* 


(?) 

Branch  points  and  cuts  of  8  and  Hq  ,  and  integration  path 
for  (12)  in  4-plane:  (a)  lossless  case,  (b)  lossy  case. 
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Fig.  5.  Ring  dipole  on  z  axis  in  infinite  half  space  (region  a) 
radiating  in  the  presence  of  contiguous  infinite  half 
space  (region  b). 


region  a  extends  from  0  to  00  and  region  b  extends  from  0  to  -00, 

while  both  regions  extend  from  to  00  in  the  x  and  y  directions. 

The  vector  potential  due  to  the  ring  dipole  is  ^-independent  and 

it  must  satisfy  the  wave  equation  in  each  region.  In  region  a  the 

vector  potential  designated  A  (p,z)  satisfies 

z 

(v~_  +  k2)  A*(p,z)  -  -u  fi(p-p')  d(z-z'),  z>0,  z ' >0  (13) 

PZ  a  Z  3. 

while  that  in  region  b  designated  Ab(p,z)  satisfies 

z 

(Vpz  +  kb)  A2(p‘z)  *  0  •  z<0  »  z’>0  (14) 

From  knowledge  of  Aa  and  Ab  the  field  components  in  the  two  regions 

z  z 

can  be  determined  in  the  usual  way  according  to 


_ / 

k2  3p3z  z 


.  co  3 
“J  u2 


(15a) 


z 


-j 


L0 

k7 

§ 


k2 

^  ft 


and 


h: 


J_  _3_ 
3p 


(15b) 


(16) 


We  observe  that  Eq.  (14)  is  a  homogeneous  differential  equation  since 
no  source  exists  in  region  b  while  the  right  hand  side  of  Eq.  (13) 
accounts  for  the  presence  of  the  ring  dipole  of  radius  p'  at  z*z'.  The 
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field  computed  from  A  and  A  must  safisfy  the  radiation  condition 

z  z 

and  the  components  which  are  parallel  to  the  interface  must  be 
continuous  at  z  =  0.  These  continuity  requirements  can  be  stated  as 


1  _  a 

k7  3p3z  z 
a 


1  _  b 

k^  3p3z  z 


at  z  =  0 


(17a) 


and 


_1_  3_  a  _  _1_  3_  b 

u  3p  z  “  y  3p  z 
a  b 


at  z  »  0  . 


(17b) 


Eqs.(13)  and  (14),  the  radiation  condition,  and  the  continuity  requirements 
(17)  completely  characterize  the  vector  potential  in  the  two  regions 
due  to  the  ring  dipole. 

We  employ  the  Hankel  transform  to  solve  (13)  and  (14)  subject  to 
(15)  and  the  radiation  condition.  Transforming  both  sides  of  (13) 
and  (14)  yields,  respectively. 


Aa(A,z)  -  -y  p'  Jn (4p ' )  6(z-z'),  z>0,  z'>0  (18) 

z  a  u 

and 

Ab(A,z)  -  0,  z<0,  z’>0  (19) 

z 

•vg  J  jj 

where  A  and  A  are  the  Hankel  transforms  of  A  and  A  and  where 
z  z  z  z 


B2  -  k2  -  A2  .  (20) 

§  § 


The  solution  of  the  differential  equation  (18)  is  obviously 
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The  first  term  of  (22)  is  the  vector  potential  due  to  the  ring  dipole 

in  medium  (y  ,e  ),  determined  as  if  the  medium  were  of  infinite 
a  a 

extent,  while  the  second  term  is  the  vector  potential  of  the  field 
resulting  from  scattering  of  the  ring  dipole  field  by  the  interface  at 
2*0. 

Turning  to  the  differential  equation  (19),  we  write  its  solution  as 
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U  jebZ 

aV,z)  =■  -j  -y  p ' Jq (Ap ' )  B(A)  ,  z<0,  z’>0 

a 


which  is  clearly  a  special  form  of  the  homogeneous  solution  of  this 

-j6bz 

equation  having  no  term  proportional  to  e  .  The  form  of  (23) 

ensures  that  the  final  solution  in  region  b  represent  a  wave  traveling 

away  from  the  interface.  Thus,  A*5  is  found  from  the  inverse  Hankel 

z 

transform  to  be 


00  1b,  z 

h  M,  .  b 

A°(P,z)  -  -j  ~  p’  A  Jq (Ap ' )  B(A)  ^ -  JQ(Ap)dA  , 


z<0 ,  z  '>  0 


To  enforce  (17b),  one  requires 


-jS  z' 

B(A)  -  C(A)  +  e  a  ,  z'>0 


which,  together  with  (24),  yields 


A^(p,z) 


00 

-j  p'  *■  J0(Ap’)  C(A) 
0 


-jB 


JQ(Ap)dA 


If,  in  addition,  (17a)  is  enforced,  one  finds  that 


-jB  z’  -jB  z’ 

8  e  3  -  C(A)  -  6,  e  3  +  C(A) 

a  b 


must  be  satisfied  from  which  C(A)  is  determined  to  be 


C  (A)  *  f  e 


-jBaz’ 


(28a) 


where 
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(28b) 


Finally,  with  C(A)  of  (28)  in  (22)  and  (26),  one  arrives  at  the 
desired  vector  potentials  in  the  two  regions: 


Aa(p,z) 

z 


00 


*  J0(*p') 


-JBjz-z'  | 
e 


JQ  (A.p)d-ft 


00  -jS  (z+z') 

-  a 

4  J0(/ip')r  - g -  JQ(*p)d*  , 

0  a 


and 

A*(p,z)  -  -j  -f  p' 

z<0,  z’>0  .  (29b) 

The  branch  points  and  branch  cuts  for  3  and  B,  ,  under  the  conditions 

a  b 

that  region  a  is  lossless  and  region  b  is  lossy,  are  illustrated  in 
Fig.  6  where  also  is  shown  the  path  of  integration  for  the  integrals 
of  (29). 

By  appealing  to  (11)  and  the  procedure  leading  to  (12),  one 
can  obtain  the  alternate  forms  for  (29)  that  appear  below. 

The  path  of  integration  for  the  integrals  of  (30)  and  (31) 
is  illustrated  in  Fig.  7. 


z>0,  z*>0 


(29a) 


~j3aZ  jB  z 

*  jQ^p'm  +  n  — g — e  J0(AP)d4  , 

a 
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Branch  points  and  branch  cuts  for  3a  and  8b  and 
integration  path  for  (29)  and  (33)  in  4-plane: 
region  a  lossless,  region  b  lossy. 


(2) 

7.  Branch  points  and  branch  cuts  for  Ba,  85,  and  Hg 
and  integration  path  for  (30),  (31),  (34),  and 
(35)  in  4-plane:  region  a  lossless,  region  b  lossy 


Green's  Function  for  Ring  Dipole  in  Region  B 


When  the  ring  dipole  of  radius  p'  is  located  at  z  =  z'  (z’<0) 
on  the  z  axis  in  region  b,  the  resulting  vector  potentials  in  the 
two  regions  satisfy 


(V*  +  kf)  A*(P*z)  "  0  »  z>0»  z  '<0 

P  Z  cl  Z 


(32a) 


and 


(VpZ  +  k^)  Ab(p,z)  =  -ub  6(p-p')  6(z-z'),  z<0 ,  z ' <0  (32b) 

plus,  of  course,  the  radiation  condition.  Also  A3  and  Ab  must  satisfy 

z  z 

the  continuity  conditions  (17).  One  may  solve  (32)  subject  to  the 
stated  conditions  by  a  procedure  paralleling  that  applied  to  (13)  and 
(14).  This  procedure  is  not  repeated  here  but  the  resulting  solutions 
are  recorded  below: 


A3(p,z)  =  -j  -j-  p' 


jV  -J8  z 

A  J0(Ap')[l-n  ~ -  e  3  JQ (Ap)dA  , 


z>0,  z'<0 


and 


(33a) 


A  z  (  P  * 2  ) 


.  Ub  , 
“J  ~2  P 


-jBJz-z' 


A  JQ (Ap ' ) 


Jq(AP )dA 


+j  -J  P' 


jB  (z+z') 
.  b 


4  j0(Ap')r  g- 


JQ (AP)dA  . 


z<0,  z'<0 


(33b) 


The  path  of  integration  of  the  above  integrals  is  the  same  as  that 
illustrated  in  Fig.  6  for  the  integrals  of  (29).  Alternate  forms  of 


48 


A®  and  are  given  on  the  next  page.  The  path  of  integration  for 
the  integrals  of  (34)  and  (35)  is  illustrated  in  Fig.  7. 
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Integral  Equations 

The  desired  integral  equations  for  the  current  on  the  ground-stake 

antenna  of  Fig.  1  now  can  be  constructed.  Due  to  the  symmetry  of  the 

structure  and  the  excitation,  all  quantities  are  independent  of  4>, 

and,  also,  the  surface  current  induced  on  the  PEC  cylindrical  tube 

is  entirely  z-directed  (J  (z)  z  ).  Under  these  conditions,  the 

— s  z  — 

induced  current  produces  an  electric  field  whose  <j>  component  is  zero 
so  the  only  component  of  electric  field  produced  by  which  is 
tangential  to  the  cylindrical  surface  of  the  antenna  is  the  z 
component.  For  convenience,  we  define  the  total  axial  currents 
Ia(z)  and  I^Cz)  on  the  portion  of  the  cylinder  in  region  a  and  region 
b,  respectively,  as 


2ttoJ  (z) 
z 


Ia(z) 

Vz) 


z  E  (0,h) 


z  E  (-d ,0) 


(36) 


where,  of  course,  a  is  the  cylinder  radius.  Furthermore,  we  define 
aa  ,  ab 

E  and  E  to  be  the  tangential  components  of  electric  field  induced 
z  z 

on  and  evaluate  at  p = a  on  the  surface  of  the  cylinder  in  region  a, 

i.e.,  ze(0,h),  due,  respectively,  to  I  and  I,,  while  Eba  and  Ebb 

a  b  z  z 

are  defined  to  be  the  z  components  of  electric  field  on  the  cylinder 
at  p  -  a  in  region  b,  i.e.,  z e  (-d,0), again  due,  respectively,  to 
I  and  I  .  Subject  to  these  definitions,  one  can  express  the  boundary 
condition  that  the  total  electric  field  tangential  to  and  on  the 
PEC  cylindrical  surface  be  zero  by  the  following  equations: 
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(37a) 


Eaa[I  ;z]  +  Eab[T;z]  =  -V  6(z-z  ), 
z  a  zb  g 


z  e  (0,h) 


[i  ;z]  +  e 

a 


bb 

z 


[Ib;z]  =  0 


z  e  (-d ,0) 


(37b) 


where  V  6  (z-z  )  accounts  for  the  delta-gap  source  of  strength  V  at 
S 

z  =  z  .  By  employing  the  vector  potential  Green's  functions  derived 
§ 

above,  one  can  write  (37)  as 


-J 


4nk 


f  >  1 

j  2 

4-t  +  k2 

dz2  a 

. 

,  > 

0 

l  0 

I  (z')  Gaa(z,z')dz'  +  [  I.(z')  Gab(z,z')dz' 
a  Jb 

-d 


*  -V  6  (z-z  ) ,  z  e  (0,h) 
8 


(38a) 


and 


-j 


4  irk. 


dz2 


+  k2 
+  Kb 


I  (z')  Gba(z,z')dz'  + 
a 


Ib(z')  Gbb(z,z')dz’ 


-d 


0, 


z  €  (-d ,0) 


(38b) 


where 


Gaa(z,z')  *  K  (z-z')  -  j 

a 


- j6  (z+z') 


4[J0(/ui)]2  r - ^ 


d*  (39a) 


(z ,  z ' )  -  -j 


j6b2’  -jS  z 

^[J0(^)]2tl-r]  — g — e  a  d* 

b 


(39b) 


Gba(z,z') 


-J  |  MJ„ 


-)6aS'  j<l. 


2  r  ,  .  n  1  e  _•*'  b 


(^)]zti+n 


e  d*. 


(39c) 
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Gbb (z , z ' )  -  K^z-z’)  +  j  |  ^[JQ(4a)]2r  - 


jp,(z+z') 

b 


d4 


(39d) 


in  which 


y~’>  ■  £ 


7T  -jk  R 
e  P 


d<f) *  ,  p*a  or  b, 


-IT 


with 


(40a) 


R2  -  (z-z')2  +  4 a2  sin2  ~  . 


(40b) 


Fundamentally,  the  left  hand  sides  of  (38a)  and  (38b)  are  special 

cases  of  (15b)  with  appropriate  Green's  functions  used  in  the  vector 

ba 

potentials.  G  and  G  are  (29a)  and  (29b),  respectively,  specialized 
to  the  present  structure  while  Gab  and  Gbb  are  specializations  of 
(33a)  and  (33b),  respectively.  The  first  term  in  each  of  (29a)  and 
(33b)  is  the  vector  potential  due  to  the  ring  dipole  in  unbounded 
space.  K&  in  (39a)  and  in  (39d)  account  for  the  first  terms, 

respectively,  of  (29a)  and  (33b).  and  of  (40)  are  employed 
here  in  (39a)  and  (39d),  rather  than  equivalent  terms  derived  from 
the  first  terms  of  (29a)  and  (33b),  simply  because  and  are 
exact  kernels  of  cylindrical  antenna  theory  and  computation  involving 
such  kernels  has  become  routine.  Of  course,  the  integrals  over  (0,°°) 
of  (39)  could  be  replaced  by  integrals  over  (-00,00)  by  appealing  to 
(30),  (31),  (34),  and  (35).  Such  a  replacement  would  allow  one 
to  convert  integrals  over  the  real  line  from  -»  to  ®  to  branch-cut 
integrals  which  in  some  cases  are  very  attractive  from  a  numerical 
computational  point  of  view. 
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4. 


COLLECTION  OF  EXPERIMENTAL  DATA 


In  this  section  is  found  a  description  of  the  apparatus  and 
procedures  to  be  used  to  obtain  data  which  will  support  the  theoretical 
analysis  of  the  foregoing  sections.  The  experimental  portion  of  this 
project  involves  two  related  efforts: 

(1)  measurement  of  the  driving  point  impedance  of  wire  antennas 
operated  near  an  air /water  interface;  and 

(2)  determination  of  the  electrical  properties  of  the  water. 
Monopole  antennas  in  various  configurations  are  to  be  positioned 

at  various  heights  above  an  air/water  interface.  The  antenna  driving 

point  is  to  be  below  the  water  surface  in  some  experiments  and  is  to 

be  above  the  water  surface,  but  behind  a  conducting  screen  oriented 

normal  to  the  water  surface,  in  other  experiments.  During  each 

experiment,  samples  of  the  water  are  to  be  gathered  and  each  sample 

is  to  be  subjected  to  several  measurement  procedures  designed  to 

determine  the  electrical  properties  of  the  water.  The  test  antenna 

configurations  and  the  techniques  for  determining  water  properties 

are  described  in  detail  later  in  this  report . 

The  experimental  effort  has  been  organized  into  three  phases: 

PHASE  I  -  Design  and  construction  of  apparatus 
PHASE  II  -  Experiments  involving  antenna  driving  point  below 
water  surface 

PHASE  III  -  Experiments  involving  antenna  driving  points  above 
water  surface,  but  behind  a  conducting  screen. 

PHASE  I  -  Design  and  Construction  of  Apparatus 

Water  Tank.  Design  and  construction  of  a  suitable  water  tank  represents 
the  most  difficult  and  expensive  task  of  this  phase.  Electromagnetic 
considerations  dictated  that  the  tank  should  be  of  non-metallic 


construction,  should  provide  a  water  surface  range  of  approximately 

2.5  m  (or  5\  in  air  at  600  Mhz),  and  should  provide  a  water  depth  of 

approximately  1.1  m  (or  20X  in  water  at  600  Mhz).  When  these 

considerations  were  balanced  against  available  space  and  materials 

and  the  fabrication  capabilities  of  this  organization,  a  general 

specification  was  made  for  a  wooden  tank  16  ft  by  16  ft  by  4  ft  -  A  in, 

overall.  Such  a  tank  would  provide  a  water  surface  4.67  m  (15  ft  - 

4  in)  by  4.67  m  and  a  water  depth  of  1.22  m  (4  ft).  (See  Fig.  1) 

Before  proceeding  with  tank  construction,  it  was  necessary  to 

verify  adequate  floor  strength  and  drainage  capacity  and  to  modify 

the  building  to  provide  a  high-capacity  water  supply  (2- in  line,  90  psi) 
» 

and  to  accomodate  suspension  of  panels  of  anechoic  materials  above 

and  around  the  water  tank.  This  was  done  by  University  of  Mississippi 

Physical  Plant  personnel  at  the  expense  of  the  School  of  Engineering. 

Details  of  tank  construction  can  be  seen  in  Figs.  1  and  2.  An 

underlying  consideration  in  the  design  of  the  tank  was  that  it  be 

capable  of  disassembly  for  economical  storage.  Accordingly,  the  floor 

was  constructed  in  four  8  ft  by  8  ft  sections  and  each  wall  was 

constructed  as  a  separate  section.  The  framework  for  floor  and  walls 

is  of  2  x  4  pine  on  16-in  centers  with  bracing  at  16-in  intervals. 

The  decking  (3/4-in  plywood)  and  the  wall  panels  (1/2-in  plywood)  were 

nailed  and  glued  to  the  framework.  The  floor  sections  were  bolted 
0 

to  the  floor  and  to  adjacent  wall  sections,  using  1/2-ii  cap  screws 
with  2-in  fender  washers  at  16-in  intervals.  Buttresses  were  added 
at  the  centerline  of  each  wall  section  to  reduce  bowing. 

In  the  original  design,  the  tank  was  to  be  made  watertight  by 
caulking  the  seams  between  pre-fabricated  sections  with  silicon  rubber 
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Fig.  2 


(a) 


.  Construction  of  floor  frame  for  water  tank. 
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1 

and  by  covering  all  wooden  surfaces  with  a  minimum  of  3  coats  of  poly¬ 
ester  resin.  In  addition,  the  edges  of  all  plywood  panels  were  coated 
with  polyester  resin  before  attachment  to  the  framework.  This  scheme 
was  not  successful,  apparently  because  the  polyester  resin  coating 
was  fractured  by  movement  within  the  structure  under  load  of  water. 

Inspection  revealed  water  penetration  into  the  plywood,  particularly 
around  nail  and  screw  heads  (even  though  the  resin  coating  had  been 
"built  up"  in  these  places)  and  along  the  manufacturer's  lamination 
seams  of  the  plywood.  The  actual  amount  of  water  lost  due  to  leakage 
was  insignificant  but  it  was  feared  that  continued  seepage  into  the 
plywood  would  cause  delamination  and  render  the  tank  unusable. 

An  attempt  was  made  to  remedy  the  seepage  problem  by  covering 
the  "trouble  spots"  with  fiberglass  cloth  submerged  in  polyester 
resin.  A  subsequent  test  revealed  that,  under  load,  water  penetrated 
the  cloth  also,  even  though  the  cloth  was  under  several  layers  of 
polyester  resin.  At  this  point,  consultation  with  the  technical 
representatives  for  a  number  of  coating  manufacturers  indicated  that 
a  plywood  tank  of  this  size  probably  could  not  be  made  entirely  water¬ 
tight  using  ordinary  coatings  (varnishes,  resins,  etc.).  The  decision 
was  then  made  to  procure  a  vinyl  liner  of  the  type  used  in  small 
swimming  pools.  A  25-mil  vinyl  liner  is  being  fabricated,  and  delivery 
is  expected  presently. 

During  the  tests  for  watertight  integrity,  it  was  determined  that 
the  tank  could  be  filled  to  full  capacity  (approximately  26000  1)  in 
less  than  one  hour.  This  corresponds  to  a  water  level  change  of 
approximately  2.5  cm  to  3  cm  per  minute. 
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Filling  is  accomplished  using  a  2-in  fire  hose  attached  to  the 
high-capacity  water  supply  (Fig.  3).  The  original  design  incorporated 
a  2-in  deck  drain  which  permitted  draining  within  l1*  hours;  however, 
installation  of  the  vinyl  liner  will  require  the  use  of  l'-s-in  drains. 

The  draining  rate  with  the  1%-in  drain  is  estimated  to  be  0.3  cm  to 
0.8  cm  water  level  per  minute.  These  rates  for  changing  water  level  are 
felt  to  be  entirely  adequate  for  conducting  the  tests  proposed. 

Slotted  Line  Apparatus.  References  [1]  and  [2]  report  favorable  results 
in  measuring  the  electrical  properties  of  water  (e  and  o)  by  means  of 
a  slotted  X-band  waveguide  apparatus.  For  this  project,  a  similar 
apparatus  has  been  constructed  from  WR137 (RG-50/U)  waveguide,  which 
is  conservatively  expected  to  provide  an  operating  band  for  dominant 
(TE^q)  mode  from  650  to  900  Mhz  when  filled  with  water.  This  apparatus 
is  depicted  in  Fig.  4.  At  650  Mhz,  the  "infinite"  section  is 
approximately  35A  in  length.  Major  components  were  obtained  from 
salvage  or  purchased  under  cost-sharing. 

In  order  to  provide  for  measurement  of  water  properties  at  lower 
frequencies  than  can  be  propagated  in  WR137  waveguide  and  to  have 
available  independently  measured  data  for  corroboration,  a  more 
elaborate  slotted  coaxial  line  has  been  designed  and  is  under 
construction.  Illustrations  of  the  design  details  and  photographs  of 
completed  parts  of  the  apparatus  are  included  in  Fig.  5.  The  coaxial 
line  is  expected  to  permit  reliable  measurements  at  frequencies  as 
low  as  500  Mhz,  at  which  the  length  of  the  "infinite"  line  is 
approximately  20A.  Introduction  of  salt  into  the  water  may  increase 
losses  along  the  "infinite"  line  to  the  point  where  frequencies  as  low 
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Fig.  3 


.  Filling  water  tank. 
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Schematic  of  slotted  waveguide  apparatus. 


Fig.  4(b).  Slotted  waveguide  apparatus  with  water  trough 
in  background. 


Fig.  4(c).  Slotted  section  in  position  for  submersion. 
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Probe  and  termination  section. 


Brass  Slide 
Contact 


67 


as  300  Mhz  may  be  used. 

For  the  TEM  propagation  mode  of  a  coaxial  line,  the  propagation 
constant  k  is  given  as 


where 


k  =  jo)/ue  =  8-ja 


e 


£-£ 

0  r 


and  u),  p,  e  and  a  have  the  usual  definitions.  To  determine  and  o 
for  the  water  dielectric,  one  proceeds  as  outlined  in  references  []] 
and  [2]  to  calculate  a  and  8  from  measured  ratios  of  the  E-field  at 
various  positions  along  the  line  to  the  E-field  at  some  reference 
position  on  the  line.  Since 


k2  -  (8-ja)2  *  (82  -  a2)  -  j2aB 


and 


then 


and 


k2  =  w2ue  =  u)2vieg£r  -  jcoyo  , 


e  _  -  -t-~  a  =  X.2  (82  -  a2) 


'r  u  ^oE0  0 


g  _  2aB  _  0  (a6v 
"  wo  “  no  17 


Each  water  sample  is  to  be  analyzed  using  both  slotted  lines.  It 
is  expected  that  the  results  will  serve  as  mutual  confirmation  of  the 
validity  of  the  techniques. 


Test  Antennas.  The  test  antennas  to  be  used  in  PHASE  II  are  sleeve-fed 
monopoles  for  which  the  radiating  structure  is  of  circular  cross-section 
and  sufficiently  small  in  radius  to  permit  thin-wire  analysis.  This 
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type  of  antenna,  when  mounted  above  a  ground  plane,  was  analyzed  by 
Butler  and  Harrison  [3]  in  1971  using  a  method  of  moments  (MOM) 
solution  to  the  Hallen's  equation.  As  part  of  this  project,  a 
MOM  solution  to  the  Pocklington ' s  equation  for  the  sleeve-fed  monopole 
has  been  developed.  This  work  was  done,  first,  to  gain  a  better 
understanding  of  these  antennas,  second,  to  obtain  baseline  data 
from  which  to  estimate  behavior  of  the  antennas  above  a  nonideal 
conducting  plane,  and,  third,  to  reflect  improvements  in  the  MQM 
technique.  Further,  since  the  performance  of  each  antenna  is  being 
verified  by  measurements  above  a  ground  plane,  the  Pocklington  code 
is  used  to  verify  antenna  performance.  Results,  to  this  point,  have 
been  entirely  consistent. 

For  each  experiment,  three  antennas  have  been  constructed,  one  of 
°  108-cm  radius,  a  second  of  0.179-cm  radius  and  a  third  of  0.318-cm 
radius.  Current  plans  are  to  conduct  some  otherwise  identical  tests 
using  each  antenna.  If  the  results  obtained  with  the  thinner  antennas 
support  the  thin-wire  assumptions  for  the  0.318-cm  antenna,  as  is 
expected,  then  the  preponderance  of  experimentation  will  be  conducted 
using  the  more  rigid  0.318-cm  antenna. 

The  sleeve-fed  monopole  antennas  have  been  fabricated  and  are  shown 
in  Fig.  6.  Note  that  the  0.318-cm  antenna  has  detachable  extensions 
to  provide  a  variety  of  antenna  heights  from  6  cm  to  25  cm  above  the 
source  gap.  The  other  two  antennas  have  fixed  heights  of  25  cm  above 
source  gap.  All  antennas  have  a  sleeve  length  of  approximately  90  cm. 
These  antennas  are  constructed  from  50-ohm  semi-rigid  coaxial  cable 
with  unjacketed  copper  outer  conductor.  The  0. 108-cm  radius  antenna 
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is  fabricated  from  0.085-in  coax,  the  0.179-cm  radius  from  0.141-in 
coax,  and  the  0.318-cm  radius  from  0.250-in  coax.  The  radiating 
structure  above  the  source  gap  is  a  hollow  tube  obtained  by  removing 
the  dielectric  and  the  inner  conductor  from  a  section  of  coax.  The 
end  nearest  the  gap  is  plugged  with  conductor,  only  to  a  depth 
sufficient  to  permit  soldering  to  the  center  conductor  from  the  feed 
sleeve.  With  ‘care  in  fabrication  and  handling,  it  is  possible  to 
have  gap  widths  of  approximately  0.5  mm,  or  about  one-half  radii  for 
the  thinnest  (0.108-cm)  antenna  to  one-sixth  radii  for  the  thickest 
(0.318-cm)  antenna.  The  driving  point  is  a  male  SMA  connector 
attached  directly  to  the  feed  sleeve. 

Straight  and  bent-wire  monopole  antennas  are  required  for 
PHASE  III;  materials  for  these  antennas  have  been  procured  and  the 
antennas  will  be  fabricated  and  tested  as  needed.  These  antennas 
are  also  to  be  constructed  from  50-ohm  semi-rigid  coaxial  cable  with 
a  male  SMA  connector  at  the  driving  point.  The  outer  conductor  will 
be  terminated  flush  with  the  front  side  of  the  conducting  screen  and 
the  center  conductor  will  be  extended  to  form  the  radiating  element. 
The  0.085-,  0.141-  and  0.250-in  coax  will  provide  radiating  elements 
with  radii  of  0.0255,  0.0455  and  0.0815  cm  respectively. 

Ground  Plane.  A  conducting  screen  is  to  be  used  during  PHASE  III 
experiments.  The  design  specifies  heavy-gauge  aluminum  to  be  clamped 
along  the  top  edge  and  suspended  as  a  curtain  along  and  just  inside 
one  wall  of  the  water  tank.  This  design  will  allow  the  aluminum  to 
be  raised  from  the  tank  when  not  in  use.  A  source  of  salvaged  1/8-in 
aluminum  panels  nominally  8  ft  by  16  ft,  but  damaged  along  the  edges. 
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has  been  located,  and  procurement  is  in  process.  Two  of  these  panels 
will  be  trimmed  and  seamed  to  form  a  ground  plane  approximately  4.5  m 
(14  ft  -  6  in)  wide  by  3  ra  (10  ft)  high. 

,  Instrumentation.  The  basic  instrument  for  measurement  of  both  the 

E-field  ratios  in  the  slotted  line  apparatus  and  the  antenna  driving 
point  impedances  is  the  Hewlett-Packard  8410A  Network  Analyzer  System. 
System  diagrams  for  both  types  of  measurements  are  shown  in  Fig.  7. 

The  measurements  process  has  been  automated  through  the  use  of  a 
Tektronix  4052  Computer/Controller .  Frequency  stability  is  maintained 
by  a  loop  containing  Hewlett-Packard  5435A  Electronic  Counter.  An 
Amplifier  Research  Model  1W1000  Broadband  Power  Amplifier  has  been 
ordered.  Incorporation  of  this  amplifier  with  its  nominal  1  W  output 
is  expected  to  boost  power  levels  adequately  to  offset  the  mismatch 
between  air-filled  and  water-filled  regions  in  the  slotted  apparatus 
and  to  accomodate  the  extended  lengths  of  antenna  feed  line  required 
when  the  driving  point  is  below  the  water  surface. 

PHASE  II  -  Experiments  Involving  Antenna  Driving  Points  Below  Water 
Surface 

.  Antenna  Configurations  and  Experiment  Geometry.  In  this  phase  of  the 

experiment,  the  antenna  configuration  is  a  straight-tube,  sleeve-fed 

•  monopole  positioned  with  its  cylindrical  axis  normal  to  the  air/water 

interface  and  with  the  source  gap  in  air  (see  Fig.  8(a)).  The 
driving  point  is  located  as  near  to  the  bottom  of  the  tank  as  is 
practical.  The  water  level  will  be  varied  to  provide  different  values 
of  gap  height  (g),  ranging  from  as  small  as  practical  to  X / 2 .  This 
experiment  will  be  repeated  for  all  three  antennas  (  a  *  0.108  cm. 
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Instrumentation  for  K-field  ratio  ■  leasurements 
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Fig.  7(b).  Instrumentation  for  antenna  impedance  measurements. 


Tank  Floor 

Fig.  8(a).  Antenna  configuration  for  PHASE  II  measurements. 


Pedestal 


Fig.  8(b).  Equipment  layout  for  PHASE  II  measurements. 
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a  =  0.179  cm  and  a  =  0.318  cm).  For  the  0.318-cm  antenna,  extenders 


will  be  used  to  provide  different  lengths  of  radiating  structure  (h) 
above  the  source  gap.  If  necessary,  these  experiments  can  be  repeated 
•for  the  thinner  antennas  by  successively  "snipping"  the  hollow  tube 
above  the  source  gap. 

Equipment  Layout  and  Procedural  Considerations.  The  antenna  is  to  be 
mounted  to  a  wooden  pedestal  such  that  the  driving  point  will  be  located 
25-30  cm  above  the  tank  floor.  This  distance  above  the  floor  is 
dictated  by  the  minimum  bend  radius  of  the  feed  line,  the  length  of 

connectors  and  adapter,  and  the  requirement  to  encase  the  connection 

« 

in  a  water-proof  seal.  The  feed  line  is  a  h-in  Spirafil  coax 
specified  to  have  losses  of  0.05  to  0.08  db/m  in  the  frequency  range 
of  interest.  Type  N  connectors  are  required  with  the  Spirafil  coax; 
therefore  an  adapter  is  required  to  accommodate  the  SMA  connector 
at  the  antenna  driving  point.  The  length  of  feed  line  will  be 
minimized  by  having  the  coax  exit  the  tank  wall  at  floor  level  and  by 
placing  the  S-pararaeter  test  set  of  the  network  analyzer  as  close  as 
possible  to  the  tank  wall.  The  calibration  technique  to  be  used  with 
the  network  analyzer  will  establish  an  electrical  reference  plane 
several  centimeters  above  the  physical  driving  point  of  the  antenna 
and  will  thus  eliminate  connector  discontinuities  from  the  final  data. 

Anechoic  materials  will  be  placed  around  and  above  the  tank  at  a 
range  of  2-3  meters  from  the  antenna.  The  antenna  will  be  located 
slightly  away  from  tank  center  to  avoid  any  degradation  of  the 
experiment  due  to  in-phase  reflections  from  the  surroundings.  Salt 
will  be  added  to  the  water  to  assure  that  the  water  is  sufficiently 


lossy  to  cause  the  currents  on  the  outside  of  the  sleeve  to  decay, 
and  hence  cause  the  sleeve  to  appear  infinite  in  length. 

A  sight  gauge  located  outside  the  tank  will  permit  accurate 
•determination  of  water  level.  In  addition,  phenomena  such  as 
capillary  action  of  the  water  on  the  walls  of  the  thin  antennas  will 
be  considered. 

PHASE  III  -  Experiments  Involving  Antenna  Driving  Points  Above  Water 
Surface,  but  Behind  a  Conducting  Screen 

Antenna  Configurations  and  Experiment  Geometry.  In  this  phase  of 
the  experiment,  the  antenna  configuration  ma'  be  a  straight-  or  bent- 
wire  monopole  positioned  above  the  air/water  interface.  The  driving 
point  will  be  located  behind  a  conducting  screen  which  intersects 
the  water  surface  at  a  right  angle  (see  Fig.  9(a)).  The  water  level 
will  be  varied  to  provide  different  separations  (h)  between  the  antenna 
and  the  water  surface.  In  some  cases  h  will  be  made  as  small  as 
practical,  which  means  that  in  experiments  with  bent-wire  antennas,  h’ 
may  be  below  the  water  surface.  As  in  PHASE  II,  some  experiments  will 
be  repeated  with  antennas  of  all  three  radii  in  order  to  validate 
the  thin-wire  assumption. 

Equipment  Layout  and  Procedural  Considerations.  The  conducting  screen 
will  be  suspended  as  near  as  practical  to  one  wall  of  the  tank  and 
lowered  until  the  bottom  edge  is  as  near  as  practical  to  the  tank 
floor.  The  antenna  mount  will  be  located  as  close  as  practical  to  the 
top  edge  of  the  tank  wall  and  offset  slightly  from  the  horizontal 
center  of  the  screen.  The  configuration  provides  screen  ranges  of 
approximately  4.5A  horizontally  and  3.5X  vertically  in  air,  and 
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Fig.  9(a).  Antenna  configurations  for  PHASE  III  measurements. 
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Fig.  9(b).  Equipment  layout  for  PHASE  I II  measurements . 
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approximately  40X  horizontally  and  20A  vertically  in  water,  at  600 
Mhz.  Greater  screen  ranges  in  air  would  be  desirable,  but  are  not 
practical  at  the  frequencies  of  interest. 

The  antennas  will  be  constructed  so  that  the  driving  point  is 
as  close  as  practical  to  the  rear  side  of  the  conducting  screen.  The 
network  analyzer  will  be  connected  directly  to  the  SMA  connector 
at  the  antenna  driving  point.  Again,  the  calibration  technique 
will  serve  to  eliminate  connector  discontinuities. 

As  in  the  PHASE  II  experiments,  anechoic  materials  will  be  placed 
around  and  above  the  tank.  Fresh  water  may  be  used  in  these  experiments 
because  increased  water  losses  are  not  required. 
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